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ABSTRACT 


Inversion  techniques  are  used  to  infer  bottom  geoacoustic  properties  using 
AN/SQQ-32  Reverberation  Level  (RL)  data  gathered  by  USS  Avenger  (MCM-1)  in 
Rhode  Island  Sound  in  February  1993.  Based  on  the  hypothesis  that  the  magnitude  of 
backscattered  energy  is  directly  related  to  the  acoustic  reflectivity  of  the  seabed,  a 
statistical  analysis  of  beam  RL  time  series  was  conducted  to  determine  its  correlation 
with  sediment  geoacoustic  character. 

A  technique  was  developed  using  the  deviation  of  the  RL  for  an  individual 
ping/beam  from  an  area-wide  average  RL  to  generate  geographic  maps  illustrating 
bottom  geoacoustic  characteristics.  Resulting  plots  of  “relative  reflectivity”  not  only 
agreed  with  the  existing  descriptions  of  sediment  distribution,  but  also  provided  more 
detailed  spatial  representation  of  bottom  geoacoustic  distribution.  This  highlighted  the 
gross  inadequacies,  particularly  in  spatial  resolution,  of  existing  information  on  bottom 
geoacoustic  distribution. 

These  plots,  when  produced  using  appropriately  small  sample  intervals,  have 
sufficient  spatial  resolution  to  expose  MCM  clutter  density  information.  Geographic 
maps  of  relative  reflectivity  can  provide  an  invaluable  aid  to  planning  search  strategy,  a 
surveying  tool  to  compare  clutter  densities,  appropriate  geoacoustic  parameters  for 
accurate  model  (including  full  wave)  predictions,  and  a  means  of  real  time  performance 
monitoring  and  assessment  (providing  the  ability  to  revise  and  modify  search  strategy). 
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I.  INTRODUCTION 


Underwater  warfare  (UW)  operations  conducted  in  shallow  water  face  a 
significant  challenge  compared  to  similar  operations  conducted  in  the  deep  open  ocean. 
The  shallow  water  littoral  regions  of  the  world’s  oceans  are  highly  complex  areas 
containing  coastal  fronts  and  eddies,  substantial  tides  and  currents,  extensive  and  often 
steep  topographic  changes,  abundant  sea  life,  varied  bottom  sediment  distributions,  and 
high  shipping  traffic.  Many  of  these  features  vary  rapidly  in  time  and/or  space  and 
combine  to  create  a  very  complicated  and  challenging  warfare  environment,  especially  to 
acoustic  sensors. 

In  shallow  water,  UW  acoustic  sensors  are  particularly  influenced  by  geoacoustic 
interactive  processes  with  the  seafloor,  which  include  transmission  into  the  sediment,  and 
reflection  and  scattering  at  the  seafloor.  Bottom  reverberation,  a  result  of  the  latter,  is 
especially  significant  and  will  often  severely  degrade  acoustic  sensor  performance  and 
capability.  As  a  result  UW  operations  are  limited  in  shallow  water  to  a  great  extent 
because  of  their  dependence  on  the  limited  information  available  on  bottom  sediment 
type  and  distribution.  Bottom  sediment  type  maps  are  based  on  large  area  averages  (low 
resolution)  and  are  typically  derived  from  a  few  in-situ  sample  locations.  Not  only  do 
they  provide  scarce  information  on  geoacoustic  parameters,  but  their  spatial  resolution  is 
usually  grossly  inadequate  for  the  high  variability  of  sediment  distributions  common  to 
shallow  water  environments. 

These  inadequacies  in  bottom  sediment  information  are  even  further  highlighted 
by  the  more  complicated  applications  involved  in  mine  countermeasures  (MCM) 
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operations.  In  addition  to  the  usual  shallow  water  limitations  faced  by  any  other  acoustic 
sensors  (reverberation,  etc.),  mine  hunting  acoustic  sensors  face  other  limitations  unique 
to  their  application.  For  example,  the  successful  detection,  classification  and 
identification  of  mines  is  extremely  dependant  on  the  number  of  false  targets,  including 
mine-like  echoes  (MILECs),  mine-like  contacts  (MELCOs),  or  non-mine  mine-like 
bottom  objects  (NOMBOs),  that  may  clutter  the  bottom  and  are  otherwise  known  as 
clutter  density  .  Also,  bottom  types  with  a  potential  for  mine  burial  may  require 
acoustic  sensors  to  have  the  capability  to  penetrate  the  sediment  in  order  to  detect  them. 
This  factor  creates  a  dilemma  for  mine  hunting,  illustrated  by  Figure  1.  While  the  higher 
frequency  acoustic  sensors  are  required  for  the  detecting,  classifying  and  resolving  detail 

of  small  targets,  lower  frequency  sensors  are  still  required  to  penetrate  the  bottom  in 
search  of  buried  targets. 

Therefore  good  high-resolution  information  on  bottom  sediment  properties  and 
distribution  is  a  major  requirement  for  all  UW  operations  in  shallow  water  and  absolutely 
essential  to  the  success  of  MCM  operations. 

The  objective  of  this  thesis  is  to  develop  a  real  time,  high-resolution  IT  method 
for  determining  bottom  sediment  acoustic  properties  and  discrete  clutter  densities  in  a 
shallow  water  environment.  This  study  sought  to  determine  those  properties  through  the 
analysis  of  an  acoustic  sensor’s  individual  beam  reverberation  level  (RL)  related  to  the 
statistical  analysis  of  many  such  RL  beam  time  series  from  the  same  area. 

The  use  of  RL  data  is  practical  not  only  because  a  large  amount  of  data  is 
a\  ailable  for  research,  but  because  the  RL  data  is  gathered  (but  currently  discarded)  by 
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active  fleet  assets  rather  than  an  oceanographic  research  vessel.  In  this  study,  the  RL  data 
was  gathered  during  an  experimental  trial,  aboard  an  Avenger  class  minehunter,  using  the 
AN/SQQ-32  advanced  minehunting  and  classification  sonar  system. 


A.  HIGH  FREQUENCY  REVERBERATION  IN  SHALLOW  WATER 

1.  MCM  Sonar 

MCM  sonars  operate  at  much  higher  frequencies,  thus  much  shorter  ranges,  than 
their  ASW  equivalents.  Most  MCM  sonars  employ  a  dual-frequency  capability; 
operating  at  lower  frequencies  (20  to  50  kHz)  to  search  and  detect  at  a  safe  stand-off 
range  (several  hundred  meters)  and  higher  frequencies  (200  to  500  kHz)  to  identify  and 
classify  targets  at  shorter  range  (Friedman,  1997;  Lathrop,  1995). 

Dedicated  US  MCM  vessels  in  use  today  are  the  Avenger  (MCM-1)  class 
minehunter  and  the  Osprey  (MHC-51)  class  coastal  minehunter,  both  equipped  with  the 
AN/SQQ-32  variable  depth  sonar  (VDS).  The  Royal  Navy’s  (RN)  MCM  vessels  are  the 
Hunt  class  minehunter/sweeper  with  the  hull  mounted  sonar  (HMS)  Type  2193  and  the 
Sandown  class  minehunters  fitted  with  the  VDS  Type  2093  (Jane’s  Underwater  Warfare 
Systems,  2000-2001). 

2.  Shallow  Water  Acoustics 

Factors  that  affect  the  transmission  of  underwater  sound  include  the  structure  of 
the  water  column,  the  interaction  at  the  boundaries  (surface  and  bottom),  and  the 
background  noise.  In  shallow  water,  it  is  the  interactive  processes  at  the  sea  bottom  that 
particularly  affect  sound  transmission,  and  often  severely.  These  processes  include 
reflection  at  the  seafloor,  transmission  into  the  sediment,  and  perhaps  the  most 
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significant,  scattering  at  the  seafloor.  Scattering  is  the  random  re-radiation  (reflection)  of 
a  portion  of  the  acoustic  energy  by  the  surface,  volume,  and/or  bottom  (Urick,  1983). 

3.  Reverberation 

Unlike  scattering,  which  is  random  in  direction,  backscatter  is  that  component 
which  is  re-radiated  (reflected)  back  towards  the  direction  that  it  came  from.  The 
combination  of  all  backscattering  mechanisms  (surface,  volume,  and  bottom)  received, 
both  directly  and  indirectly,  by  the  acoustic  monostatic  (co-located  transmitter/receiver) 
sensor  of  origin  is  known  as  reverberation.  Although  subtly  different  in  their  definition, 
backscattering  and  reverberation  are  directly  related  and  often  considered  the  same. 

Sea  surface  reverberation  is  due  to  backscattering  from  a  roughened  air-ocean 
interface,  a  function  of  wind-driven  wave  heights,  and  therefore  directly  related  to  wind 
speed.  Volume  reverberation  occurs  due  to  the  backscattering  from  marine  life  or  other 
suspended  material  within  the  water  column  and  from  inhomogeneities  in  the  ocean 
thermal/density  structure.  Bottom  reverberation  is  typically  the  most  dominant  and 
variable  component  of  reverberation  in  shallow  water,  and  is  due  to  the  backscattering 
from  the  ocean-bottom  interface. 

In  shallow  water,  the  effects  of  the  bottom  become  particularly  significant  due  to 
the  nearness  of  the  seabed.  The  reflective  nature  of  the  overlying  sediment,  and 
sometimes  the  underlying  material,  can  exert  a  significant  effect  on  active  sonar 
performance.  Resulting  reverberation  levels  (RL)  can  be  overwhelming  to  the  point  that 
a  target’s  echo  level  (EL)  is  saturated  or  masked  by  the  RL,  and  hence  not  detected. 
Acoustic  sensor  operations  in  the  shallow  water  environment  are  further  complicated  by 
the  high  spatial  variability  of  sea  bottom  composition.  Clearly,  a  detailed  knowledge  of 
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the  bottom  environment  in  shallow  water  and  how  it  effects  bottom  reverberation  is 
critically  important  and  essential  to  detection  ranges  and  their  accurate  prediction. 
Bottom  backscattering  processes,  and  hence  bottom  RL,  are  complexly  inter-dependent 
on  frequency,  grazing  angle,  and  the  nature  of  the  bottom  (both  grain  size  and  surface 
roughness). 

a)  Backscattering  Strength  (Sb) 

The  intensity  of  the  backscattered  energy  from  the  sea  floor  that  arrives  at 
the  sonar  receiver  is  parameterised  as  the  bottom  backscattering  strength  (Sb)  (Urick, 
1983,  pp.  238-239); 

Sb=101ogi-  , 

Ainc 

where  Isca,  is  the  intensity  of  the  scattered  sound  per  unit  area  at  a  unit  distance  and  Iinc  is 
the  incident  plane  wave  intensity  impinging  on  the  scattering  area. 

b)  Grazing  Angle  ( 6) 

Grazing  angle  is  the  horizontal  arrival  angle  measured  from  the  seabed 
(whereas  the  incident  angle  is  measured  from  normal  to  the  seabed).  Figure  2  illustrates 
the  scattering  strength  (Sb)  as  a  function  of  grazing  angle.  It  indicates  that  scattering 
strength  tends  to  be  independent  of  grazing  angle  at  small  angles,  with  increasing  slope 
from  8°  to  20°  where  Bragg-type  scattering  occurs,  and  then  rising  steeply  as  the 
scattering  approaches  coherent  reflection  from  the  seafloor. 

c)  Frequency 

Backscattering  coefficients  are  independent  of  frequency  when  the  spatial 
scale  of  bottom  roughness  is  large  compared  to  the  wavelength.  However  when  the  scale 
of  bottom  roughness  (appreciable  roughness  at  small  spatial  scales)  is  small  compared  to 
the  wavelength,  the  scattering  strength  increases  with  frequency.  Thus  variation  of 
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bottom  backscattering  with  frequency  is  highly  complex  because  it  is  influenced  by  both 
bottom  composition  and  bottom  roughness  (AOMC  Training  Publication). 

d)  Grazing  Angle  ( 6)  and  Bottom  Type  ( Roughness ) 

Measurements  of  RL  and  Sb  arising  from  seabed  backscattering  have 

indicated  that  Sb  is  a  function  of  grazing  angle  (0)  and  bottom  type.  Backscattering 

strength  from  acoustically  rough  topography  is  traditionally  represented  using  a 
Lambert’s  law  relationship; 

Sb  =  10  log  p.  +  10  log  sin20, 

where  p,  Mackenzie’s  constant,  represents  the  degree  of  reflectivity  of  the  sediment 
interface.  Based  on  the  deep  water  data  of  Mackenzie  (1961)  and  others,  Urick  (1983, 
pp.  27 1  -273)  reported  that  over  moderately  rough  topography  of  varying  sediments  and 
small  to  intermediate  grazing  angles,  Mackenzie’s  constant  (k  =  10  log  p)  averaged  about 
-_9  dB.  and  was  essentially  invariant  with  frequency  and  grazing  angle.  More  recent 
shallow  water  experiments  (south  of  Long  Island)  have  indicated  an  association  of 
Mackenzie's  constant  to  the  bottom  sediment  type  and  that  grazing  angle  dependence 
ma.v  be  scaled  as  sin  0  over  muddy,  absorptive  bottoms  where  volume  reverberation  from 
the  upper  layers  of  the  sediment  may  be  important  (Scanlon  et  al,  1995;  McCammon, 
1993).  Exercise  LWAD  99-1  reproduced  these  relationships,  with  measured  bottom 
backscattering  strengths  over  a  silty  clay  bottom  having  a  sin  0  relationship  and  that  over 
a  sandy  bottom  exhibiting  a  sin20  relationship  (Schalm,  1999).  These  findings  suggested 
a  direct  correlation  exists  between  backscattering  strength  and  sediment  type,  and  with 
this  knowledge,  Schalm  (1999)  used  inverse  techniques  to  infer  the  nature  of  the 
sediment  type  as  a  function  of  its  reflectivity  from  RL  data. 
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e)  Grain  Size 

The  surface  sediment  type  can  also  be  characterised  by  its  grain  size. 
Various  texts  and  papers  outline  such  classifications,  all  having  similar  categorisations, 
for  example  Komar  (1998).  Although  a  relationship  between  grain  size  and  Sb  of  graded 
sands  has  been  demonstrated  in  the  laboratory  (Nolle,  1963),  a  corresponding  relationship 
has  not  been  verified  in  the  field  (Figure  3).  Statistical  examination  of  most  available 
data  suggests  that  there  are  at  least  two  independent  scattering  processes;  a  scattering 
process  that  is  driven  by  grain  size  and  another  mainly  the  result  of  scattering  from  gas 
bubbles  trapped  within  the  sediment  (Chotiros,  1992). 

B.  GEOACOUSTIC  PREDICTION 
1.  Traditional  Approach 

Traditionally  in  deep  water  where  bottom  sediment  properties  are  relatively 
constant  spatially,  bottom  backscatter  coefficients  have  been  determined  from  measured 
data  at  single  locations  (Gomes  et  Ah,  1998;  Urick,  1983;  Boehme  et  ah,  1988;  Boehme 
et  al,  1985;  McKinney  and  Anderson,  1964;  Jackson  et  ah,  1986;  Stanic  et  ah,  1988; 
Stanic  et  ah,  1989;  Urick,  1954;  Hamilton,  1980;  Turgut  et  ah,  1990;  Nolle,  1963; 
Medwin,  1977).  Estimates  of  bottom  backscatter  coefficients,  derived  from  these  spot 
measurements,  have  then  been  applied  over  large  geoacoustic  provinces  in  “look  up 
tables”  for  application  in  deep  water  acoustic  propagation  models. 

In  shallow  water  this  application  is  particularly  inappropriate  as  the  bottom 
sediment  types  and  clutter  densities  are  extremely  variable  on  spatial  scales  of  hundreds 
of  meters,  not  hundreds  of  miles.  Scanlon  et  al.  (1996)  showed  that  the  spatial  scale  can 
grade  from  hard  sand  to  mud/silt  clay  sediments  over  distances  less  than  500  m.  This  has 
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been  supported  by  numerous  very  high  frequency  side  scan  sonar  (SSS)  measurements 
using  the  AQS-14  or  Klein  5000  SSS  (Null  et  al.,  2000).  Thus  averaging  sonar  data  over 
large  areas  in  shallow  water  to  determine  an  average  bottom  backscatter  coefficient 

produces  gross  inaccuracies  in  bottom  effects  and  hence  unreliable  predictions  of  sonar 
performance. 

Another  limitation  to  this  approach,  of  using  backscatter  coefficients,  is  the 
assumption  that,  at  the  high  frequencies  associated  with  MCM  sonars,  the  scattering 
occurs  at  the  sediment-water  interface.  Recent  studies  (Neumann  et  al.,  1999)  show  that 
the  scattering  occurs  from  inhomogeneities  in  the  sediment  sub-bottom  for  slow  (i.e., 
clay,  mud)  sediments  at  low  grazing  angles. 

2.  Alternative  Approach 

Although  widely  spaced,  single-point  bottom  backscattering  measurements  cannot 
accommodate  or  resolve  the  high  spatial  sediment  variability  inherent  in  shallow  water, 
acoustic  inverse  techniques  (IT)  can  be  developed  to  estimate  this  variability,  based  upon 
an  analysis  of  MCM  sonar  bottom  reverberation  observations.  The  IT  is  first  used  to 
determine  the  type  and  spatial  distribution  of  sediments  on  small  spatial  scales,  as  well  as 
discrete  clutter  (MDLECs,  MILCOs,  and/or  NOMBOs).  A  bottom  loss  model  can  then  be 
matched  to  the  highly  variable  bottom  in  the  shallow  water  area  “surveyed”  with  the  IT 
estimate  derived  from  the  acoustic  data. 

If  SSS  bottom  mapping  data  are  available,  then  the  lower  frequency  search  sonar 
IT  estimates  can  be  fused  with  the  SSS  data.  High  frequency  SSS  data  provide  a  very 
high  spatial  resolution  of  bottom  surface  features,  but  contain  little  information  regarding 
buried  features  such  as  mines  and/or  buried  clutter  (MELECs,  MILCOs,  and/or 
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NOMBOs),  not  to  mention  the  vertical  gradient  of  sediment  properties.  As  Figure  1 
shows  schematically,  lower  frequency  MCM  and  ASW  sonars  can  provide  sediment- 
penetrating  capabilities,  but  have  less  seafloor  surface  resolution  than  higher  frequency 
SSS  data.  Although  multi-sensor  data  fusion  capability  is  still  in  development, 
“surveying”  the  same  bottom  area  with  more  than  one  sonar  is  highly  desirable. 

Once  the  IT-derived  sediment  type  and  clutter  objects  are  mapped  (and  archived), 
one  may  perform  propagation  and  sonar  performance  modelling  in  one  of  two  ways.  The 
classical  bottom  backscatter  “look  up  tables”  can  be  input  into  the  propagation  model 
based  on  IT  results  or  one  can  use  a  full  wave  propagation  model  without  look  up  tables. 

t 

C.  INVERSION  TECHNIQUES  (IT) 

In  ocean  acoustics  the  "forward  problem"  addresses  those  methods  which  predict 
the  acoustic  propagation  characteristics  based  on  a  knowledge  of  the  relevant  oceanic 
parameters  (e.g.,  sound  speed  structure,  sea  surface  roughness,  geoacoustic  properties, 
etc.).  An  "Inverse"  (or  inversion)  technique  is  the  exact  opposite  such  that  one  infers  one 
or  more  oceanic  parameters  through  the  analysis  of  acoustic  propagation  data.  (Medwin 
etal.,  1998). 

Inverse  techniques  have  a  wide  variety  of  applications  including  the  simulation  of 
broadband  signal  propagation  (Westwood  et  al.,  1987)  and  high-resolution  geoacoustic 
modelling  (Osier  et  ah,  1997).  ITs  have  also  been  used  to  identify  features  in 
reverberation  data  from  which  acoustic  properties  such  as  sound  speed  structure,  surface 
roughness,  and  geoacoustic  properties  can  be  determined.  Schalm  (1999)  used  such  an 


9 


IT  to  infer  bottom  sediment  characteristics  from  analysing  beam  reverberation  level  (RL) 
data  by  an  active  ASW  sonar  (AN/SQS-53C). 

A  similar  IT  is  used  in  this  study  to  identify  features  in  beam  reverberation  level 
(RL)  time  series  of  active  MCM  sonar  (AN/SQQ-32)  data  from  which  the  bottom 
acoustic  properties  are  deduced.  Specifically,  features  of  the  RL  time  series  and 
statistical  analysis  describing  the  character  of  the  RL  time  series  are  analysed  to  identify  a 
correlation  to  sediment  type. 

As  previously  discussed  (Chapter  I  B.2),  reverberation  is  a  combination  of 
backscattered  energy  from  the  sea  surface,  within  the  water  column  (volume),  and  the  sea 
floor.  For  this  study  it  is  assumed  that  the  volume  and  sea  surface  components  of 
reverberation  were  constant  within  the  temporal  and  spatial  scales  involved.  Therefore,  it 
is  hypothesised  that  for  a  single  ping  from  the  AN/SQQ-32  sonar  ,  the  individual  RL 
beam  (bearing)  time  series  should  only  differ,  from  beam  to  beam,  because  of  variations 
in  backscattered  energy  from  different  bottom  types.  This  hypothesis  is  further  extended 
to  encompass  a  number  of  pings  (as  many  as  442  pings  were  used  in  one  combined  data 
set >.  as  long  as  they  occur  in  the  same  general  area,  and  that  any  changes  to  the  ping  type 
or  mode  of  operation  (source  level,  pulse  length,  etc.)  are  accounted  for. 

The  following  assumptions  were  made  to  support  this  hypothesis: 

(1)  The  sound  speed  profile  in  the  area  around  the  source  (within  ~3  km)  was 
constant  over  the  period  of  observations  (12  minutes). 

(2)  For  each  ping  transmission  the  output  source  level  (SL)  from  the  AN/SQQ-32 
in  any  direction  did  not  vary  (i.e.,  the  SL  was  independent  in  azimuth). 
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(3)  There  were  no  major  bathymetric  or  topographic  changes  that  would 
significantly  change  the  RL  curve  (i.e.,  no  upslope/downslope  effects). 
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II.  DATA 


A.  DATA  COLLECTION 

The  full  resolution  reverberation  level  (RL)  data  examined  in  this  study  was 
collected  by  the  Applied  Research  Laboratory  -  University  of  Texas  (ARL-UT)  aboard 
USS  Avenger  (MCM-1)  on  19  February  1993  in  Rhode  Island  Sound  over  the  New 
England  continental  shelf  (Figure  4).  The  Avenger  operated  within  an  exercise  minefield 
approximately  8  km  SE  of  Newport,  RI,  in  25-27  m  depth  of  water  (Figure  5).  The 
minefield  consisted  of  15  exercise  mines  of  varying  types,  both  bottom  and  moored 
(Figure  6).  Over  600  pings  from  USS  Avenger’s  AN/SQQ-32  Mine-Hunting  sonar  were 
recorded  between  1 1:02  and  13:07  (Eastern  Time)  over  two  identical  navigational  tracks 
(course  approximately  270°),  within  the  vicinity  of  the  minefield  (Figure  6).  These  pings 
were  transmitted  along  a  distance  of  1500  m  on  Track  1,  and  900  m  on  Track  2.  The 
orientation  was  such  that  the  ping  locations  from  the  last  400  m  of  Track  1  overlapped  the 
ping  locations  of  the  first  400  m  of  Track  2  (Figure  7). 

1.  Weather  and  Oceanography 

All  acoustically  relevant  meteorological  and  oceanographic  parameters  were 
reported  as  relatively  invariant  for  the  exercise,  making  for  near-laboratory  conditions. 
The  wind  speed  was  less  than  5  ms'1  from  the  NW,  the  sea  surface  was  smooth  with  a 
significant  wave  height  of  less  than  0.3  m  (Pitt,  2000).  However,  meteorological 
observations  from  two  buoys  in  the  vicinity  of  Rhode  Island  Sound  (Table  1)  showed  that 
the  mean  wind  speed  gradually  decreased  from  4-6  ms'1  to  3-4  ms'1  over  the  two  hour 
duration  of  the  exercise  (NDBC,  2000).  Such  a  reduction  is  unlikely  to  have  affected 
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acoustic  conditions  within  the  duration  of  each  track,  but  the  decrease  could  have  been 
sufficient  to  make  a  significant  difference  to  the  sea  state  (surface  roughness)  experienced 
between  tracks  one  and  two  (separated  by  over  1.5  hours). 


Buoy  Location 

Buzzard’s  Bay 

Long  Island  Sound  1 

Local  Time 
2/19/1993 

Wind 

Direction 

Wind  Speed 
Mean/Gusts  (ms'1) 

Wind 

Direction 

Wind  Speed 
Mean/Gusts  (ms'1) 

11:00 

330 

5.8/ 7.5 

332 

4.3/ 5.6 

12:00 

334 

4.9/ 5.4 

318 

4.1/ 5.6 

13:00 

332 

3.7/ 4.2 

311 

4.0/ 4.9 

Meteorological  Data  from  Buoys  in  the  Vicinity  of  Rhode  Island  Sound. 


Climatological  data  suggests  that  the  mixed  layer  extended  to  the  seabed  and  this 
was  confirmed  by  observation  taken  on  19  February  (Pitt,  2000).  The  temperature  profile 
was  essentially  isothermal  representing  typical  oceanic-atmospheric  forcing  for  these 
waters  during  winter.  Little  spatial  or  temporal  deviation  is  expected  within  the  confines 
of  the  study  area.  The  resulting  positive  sound  speed  profile  (SSP)  provides  a  slight 
upward  refractive  effect  to  sonar  transmissions,  which  can  be  considered  negligible 
within  the  depths  and  ranges  of  this  study.  The  wave-guide  is  so  shallow  (25-27  m)  that 
the  declination  angle  rays  essentially  interact  with  the  boundaries,  reflecting  from  the  sea 
surface  and  ocean  bottom.  Direct  paths  within  the  water  column  at  declination  angles 
existed  only  for  very  short  ranges  and  large  incident  angles.  Their  grazing  angles  ranged 
from  approximately  1°  at  1100  m  range  to  approximately  15°  at  80  m  range,  all  defined 


as  shallow  angles  from  Chapter  1 . 
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2. 


Sediments  and  Geoacoustics  of  Area 


a)  The  New  England  Continental  Shelf 

Most  of  the  continental  shelf  south  of  New  England  is  covered  with  sands 
and  locally  gravelly  sands  (Figure  8).  The  finer-grained  materials  in  these  sediments 
have  been  winnowed  out  and  transported  either  shoreward  into  estuaries  or  off  the  shelf 
and  into  canyons  onto  the  continental  slope  (Tucholke,  1987).  Studies  by  Milliman 
et  al.,  (1972)  found  most  of  the  coarse  fraction  (125-250  mm)  is  comprised  of  quartz  and 
feldspar,  although  there  are  local  concentrations  of  glauconite  and  phosphorite.  Rivers 
draining  from  the  eastern  part  of  North  America  carry  little  sediment  to  the  continental 
shelf  (Milliman  et  al.,  1983);  most  of  that  sediment  is  trapped  in  estuaries  or  coastal 
marshes  (Tucholke,  1987).  Relict  sediments,  those  left  from  a  glacial  deposition  cycle, 
are  concentrated  on  the  continental  shelf  and  consist  of  fluvial  and  carbonate  detritus  and 
glacial  outwash. 

In  summary  shelf  sediments  of  this  region  are  generally  comprised  of 
mixed  subarkosic  to  arkosic  sands  of  glacial  origin  with  10-25%  feldspar  /  feldspar  and 
quartz  (Tucholke,  1987). 

b)  Rhode  Island  Sound  and  the  Area  of  Study 
An  extensive  survey  of  Rhode  Island  Sound  surface  sediments  by  McMaster 

(1960)  revealed  that  the  region,  in  which  the  area  of  this  study  is  located,  consists  of  well- 

sorted  fine  sand.  A  large  area  of  silty  sand  occurs  adjacent  to  the  Rhode  Island  mainland 

and  a  tongue  of  this  finer  sediment  lies  approximately  3  km  WNW  of  the  study  area.  The 

MIW  pilot  (NAVO)  described  the  area  as  “...covered  by  irregular  interwoven  patches  of 

gravel,  sand  and  mud.  It  is  often  strewn  with  rocky  outcrops  that  sometimes  contain 
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boulders.  Seaward  of  the  sounds,  the  bottom  is  smoother  and  consists  of  sand  and 
sand/gravel  mixes.” 

Spatially  detailed  sediment  information  for  the  specific  area  of  study  was 
unavailable.  However,  two  sources  did  provide  some  indication  of  the  spatial  variability. 
The  ‘Bottom  Materials’  chart  from  the  NAVO  MIW  Pilot  (2000)  of  this  region  (Figure  9) 
implies  that  most  of  the  area  is  comprised  of  mixed  sand,  gravel  and  rocks,  but  may 
include  mud  and  sand  in  the  northeast  part  of  the  study  area.  It  also  suggests  that  the  far 
west  of  the  study  area  is  more  likely  to  be  just  sand.  The  ‘nature  of  the  bottom’ 
descriptors,  detailed  on  the  navigational  chart  (#13221)  of  this  area  (Figure  10),  generally 
supports  this  picture  with  sands  throughout,  rocks  and  pebbles  in  the  eastern  half  of  the 
area,  and  an  adjacent  area  to  the  north  is  described  as  ‘hard’.  In  addition,  high-resolution 
bathymetry  obtained  for  this  area  (Figure  11)  shows  little  depth  variation  (<  1  m) 
implying  that  the  bottom  geomorphology  is  relatively  featureless  especially  in  the  west. 
However,  an  area  in  the  middle  appears  relatively  more  erratic  and  could  certainly  be 
accounted  for  by  rocks  or  rock  outcrops  implied  in  the  descriptors  above. 

Unfortunately,  there  were  no  sediment  cores  available  for  this  particular 
area,  w ith  the  nearest  such  cores  all  over  5  km  away. 

3.  Description  of  AN/SQQ-32  Reverberation  Data 

The  beam-formed  sonar  reverberation  level  data  used  in  this  study  was  recorded 
in  full  resolution  from  USS  Avenger’s  Raytheon  /  Thomson  Sintra  AN/SQQ-32  sonar. 
Operation  was  in  an  experimental  search  mode  at  a  frequency  range  of  35  to  36  kHz.  The 
search  mode  has  a  90°  azimuthal  transmit  coverage  (Pitt,  2000).  The  receiver  employs  28 
beams,  each  with  a  2.5°  beamwidth  giving  an  azimuthal  receive  coverage  of  70°.  This 


16 


azimuthal  beam  pattern  is  duplicated  vertically  to  produce  two  different  azimuthal 
coverage  sectors  at  upper  and  lower  depression  angles.  The  depression  angles  throughout 
this  exercise  were  set  at  4°  and  12°,  respectively  (Pitt,  2000).  Data  from  the  lower 
depression  angle  of  -12°  was  used  for  analysis  in  this  study  to  ensure  strong  bottom 
interaction  of  the  acoustic  signal  despite  the  slight  upward  refractive  effect  of  the  water 
column. 

During  the  exercise  the  SQQ-32  sonar  used  all  three  of  its  transmission  power 
modes:  low,  medium  and  high.  The  low  and  medium  power  transmissions,  210  dB  and 
215  dB  source  levels,  respectively  (peak  re  pPa  @  1  m),  had  an  experimental  pulselength 
of  11  ms  LFM,  while  the  high  power  transmission  (220  dB)  had  an  experimental 
pulselength  of  2  ms.  All  three  modes  of  operation  had  a  bandwidth  of  9  kHz.  Pings  were 
transmitted  at  a  rate  of  approximately  one  every  2  seconds  along  Track  1  and  one  every 
second  along  Track  2  (Pitt,  2000). 


Track 

Event 

Start  Time 

End 

Time 

Number 
of  Pings 

Pulse 

Length 

Power 

Setting 

Source 

Level 

1 

1 

109 

1 1ms 

Medium 

215dB 

1 

2 

11:05:26 

11:12:20 

1 1ms 

Low 

210dB 

1 

3 

88 

1 1ms 

Medium 

215dB 

4 

12:45:09 

12:49:36 

156 

2ms 

High 

220dB 

Table  2.  Details  of  Recorded  AN/SQQ-32  Transmissions  on  19  February  1993. 

Each  period  for  which  pings  were  transmitted  at  a  particular  power  setting  is 
referred  to  as  an  event.  These  events  are  detailed  in  Table  2.  Events  1  through  3 
occurred  on  Track  1;  Event  1  (Medium  Power)  included  109  pings  covering  the  first  400 
m.  Event  2  (Low  Power)  included  228  pings  covering  the  central  800  m,  and  Event  3 
(Medium  Power)  included  88  pings  covering  the  last  300  m  (Figure  12).  Event  4  (High 


17 


Power)  included  156  pings  and  exclusively  occurred  along  the  entire  900  m  of  Track  2. 
Throughout  all  events  the  sonar  heading  was  the  same  as  ship’s  head  (~270°T),  with  the 
exception  of  some  45  pings  at  the  end  of  Event  4  which  progressively  veered  from 
~270°T  to  ~300°T  as  the  sonar  was  trained  to  starboard.  The  approximate  range  of 
coverage  for  the  pings  analysed  in  this  study  was  1 100  m  giving  a  total  study  area 
approximately  3000  m  long  and  1000  m  wide. 

4.  Description  and  Preparation  of  Data 

The  beam-formed  sonar  reverberation  data  from  all  601  pings  was  originally 
recorded  at  a  sample  rate  of  18  kHz  and  supplied  in  binary  format,  along  with  ASCII  log 
files,  all  on  9  CD  ROMs.  Positional  information  was  not  included,  but  was  obtained  by 
matching  ping  times  to  the  ship’s  GPS  navigation-track  log  output  supplied  separately  in 
ASCII  format.  The  data  was  first  converted  from  binary  to  ASCII  files  (‘.mat’)  for 
analysis  using  Matlab®.  For  each  ping  the  data  was  then  processed  to  convert  from  a 
complex  measure  of  intensity  (I)  in  volts  (v)  to  reverberation  level  (RL)  in  decibels  (dB 
re  1  pPa  @  lm).  This  was  achieved  using  the  following  equation  (Pitt,  2000): 

RL  =  10  x  logl0(Abs(I))  -  TVG  -  RS 

where  TVG  =  time  varying  gain,  and  RS  =  receiver  sensitivity  (-69dB  re  1  pPa  @  lm). 
The  Time  Varying  Gain  (TVG),  applied  at  the  preamps  before  recording,  is  detailed  in 
Table  3  and  illustrated  in  Figure  13.  It  was  removed  (subtracted)  from  the  data  in  the 
same  way  it  was  applied. 
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Range  (yds) 

Gain  (dB) 

Range  (yds) 

Gain  (dB) 

Range  (yds) 

Gain  (dB) 

85-95 

-5 

280-325 

13.0 

705  -  750 

28.0 

95-115 

-2 

325  -  360 

14.5 

750-805 

29.5 

115-130 

-0.5 

360  -  390 

31.0 

130-150 

1.0 

390  -  425 

17.5 

32.5 

150-160 

2.5 

425  -  455 

19.0 

915-965 

34.0 

160  -  195 

5.5 

455  -  490 

20.5 

965  -  985 

35.5 

195-215 

7.0 

490  -  575 

22.0 

985  -  1030 

35.5 

215  -  230 

8.5 

575  -  605 

23.5 

1030  -  1095 

37.0 

230  -  255 

10.0 

605  -  655 

25.0 

1095-  1165 

38.5 

255  -  280 

11.5 

655-705 

26.5 

1165  -  end 

40.0 

Table  3.  Gain  vs  Range  of  TVG  Function.  (From  Pitt,  2000) 

For  each  ping  the  processed  data  consisted  of  56  time  series,  one  for  each  of  the 
28  beams  at  both  upper  and  lower  declination  angles.  For  each  beam  the  time  series 
contained  approximately  25,000  sample  points,  which  translates  to  about  1.5  seconds  in 
time  or  1 100  m  in  range  (assuming  two-way  travel  at  1500  ms'1  sound  speed).  Figure  14 
displays  a  typical  processed  RL  time  series  from  a  single  beam  (beam  14).  The 
significant  change  in  slope  of  RL,  which  occurs  at  approximately  0.2  seconds,  represents 
the  transition  from  spherical  to  cylindrical  spreading  due  to  boundary  interaction.  The 
relatively  unchanging  slope  of  RL  beyond  0.2  seconds  implies  that  the  RL  decays  with 
range  at  a  fairly  constant  rate  due  to  cylindrical  spreading,  suggesting  that  there  is  little 
change  in  the  environment  over  this  range. 


B.  PROCESSING  METHOD 

An  example  of  how  the  sediment  effects  reverberation  can  be  seen  by  comparing 

two  different  RL  time  series  in  which  distinctly  different  characteristics  can  be  seen.  The 

single  beam  (beam  14)  RL  time  series  in  Figure  14  shows  that  the  RL  curve,  beyond  0.2 
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seconds  decays  at  a  fairly  constant  rate  due  to  cylindrical  spreading.  In  contrast,  Figure 
15,  another  single  beam  (beam  21)  RL  time  series  from  a  different  ping,  shows 
significant  changes  in  the  slope  of  RL  at  0.6  and  0.75  seconds,  well  beyond  the  spherical 
to  cylindrical  spreading  transition  range  (~0.2  seconds).  In  fact,  it  is  clear  that  the  RL 
actually  increased  during  this  period.  This  would  imply  that  the  environment  changed 
significantly  and  this  heightened  degree  of  backscatter  is  due  to  the  presence  of  a  more 
reflective  bottom. 

The  inversion  technique  (IT)  proposed  for  RL  curve  analysis  is  as  follows:  (1) 
determine  a  characteristic  RL  curve  for  an  area  determined  from  the  average  of  many 

pings  and  (2)  compare  individual  beam  data  to  determine  its  difference  from  the  area¬ 
wide  mean. 

A  variety  of  statistical  parameters  derived  from  the  beam  RL  time  series  can  be 
examined  to  determine  if  the  characteristics  of  the  RL  curve  bear  any  correlation  to 
sediment  type.  Schalm  (1999)  investigated  the  mean  slope,  standard  deviation,  skew,  and 
excess  kurtosis,  but  only  was  only  able  to  show  promising  results  with  the  mean  slope 
due  to  time  and  funding  limitations.  Therefore  the  mean  RL  and  the  slope  of  RL  were 
investigated  here  to  build  on  Schalm’s  work. 

Before  continuing,  some  background  about  the  statistical  method  used  in  this 
stud\  is  required.  The  statistics  calculated  are  not  based  on  the  entire  time  series  but 
rather  upon  segments  of  the  time  series  partitioned  into  intervals  of  50,  100,  200,  500,  or 
1000  sample  points  corresponding  to  2,  4,  8,  20,  or  40  meters  range  spacing,  respectively. 

In  addition,  a  degree  of  overlap  of  the  data  is  often  useful,  e.g.,  25,  50  or  75%  of  the 
interval  and,  based  on  the  results  of  Schlam  (1999),  a  75%  overlap  was  used.  The 
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statistics  derived  from  the  overlapped  segments  are  referred  to  as  "moving"  statistics 
(moving  average  or  mean).  For  example,  a  1.4-second  time  series,  with  18,000  data 
points  per  second,  contains  25,200  data  points.  With  an  interval/overlap  combination  of 
1000/75%  the  first  data  point  of  the  moving  mean  would  be  the  mean  of  points  1  to  1000, 
the  second  point  would  be  the  mean  of  points  251  to  1250,  the  third  would  be  the  mean  of 
points  501  to  1500,  etc.  Selection  of  interval  entirely  depends  on  the  degree  of 
smoothing  desired  and  the  spatial  size  of  the  features  of  interest.  Based  mainly  on  limited 
data  capacities  an  interval  of  1000  was  selected  and  used  as  the  standard  throughout  the 
analysis.  Figure  16  shows  an  individual  single  beam  (beam  14)  RL  time  series  (black) 
overlaid  by  the  1000/75%  moving  mean  (white)  RL  curve. 

1.  Example  Using  RL  Data  at  1000/75%  Moving  Mean 

From  the  88  pings  of  Event  3  an  ‘Event-wide’  average  RL  curve  is  calculated  by 
first  taking  the  1000/75%  moving  mean  for  all  of  the  beams  (28)  in  each  transmitted  ping 
to  smooth  the  curve.  (The  slope  of  each  beam  can  be  determined  by  calculating  the 
gradient  (d/dt)  of  the  smoothed  RL  curve).  These  smoothed  curves  from  each  beam  are 
then  used  to  calculate  the  mean  RL  curve  for  that  particular  ping.  This  process  is 
repeated  for  all  the  pings  in  Event  1  generating  88  individual  mean  RL  curves,  each 
representing  1  ping  over  the  70°  sector  width.  The  ‘Event-wide’  RL  curve  is  then 
calculated  from  all  the  pings’  individual  mean  RL  curves.  Figure  17  shows  all  the  ping’s 
individual  mean  RL  curves  for  Event  3  (black)  overlaid  by  the  Event-wide  mean  RL 
curve  (white).  Calculating  the  mean  from  all  (24x88)  RL  data  sets  from  Event  3 
generates  an  Event-wide  average  RL  curve  based  on  over  2000  individual  RL  time  series. 
This  process  was  performed  on  the  data  sets  of  each  of  the  four  separate  events. 
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Figure  18  shows  the  smoothed  RL  curve  for  a  single  individual  beam  (beam  14) 
(dashed  line)  along  with  its  appropriate  Event-wide  average  RL  curve  (solid  line).  Each 
point  along  the  smoothed  single  beam  RL  curve  has  an  associated  range  and  bearing 
which  corresponds  to  a  geographic  position.  By  determining  the  difference  between  the 
area-wide  average  RL  curve  and  the  individual  smoothed  beam  RL  curve,  a  geographic 
contour  plot  can  be  generated  to  display  the  anomaly  from  the  Event-wide  mean  (Figure 
19).  Subsequently,  different  pings  can  be  compared  to  determine  if  any  geographic 
correlation  exists  between  sets  of  pings  to  identify  a  unique  but  consistent  area  exhibiting 
similar  RL  characteristics.  In  Figure  19,  red  and  yellow  represents  areas  of  relatively 
high  reflectivity  while  green  and  blue  are  relatively  low  reflectivity  areas. 
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III.  ANALYSIS 


A.  ANALYSIS  OF  RL  DATA 

1.  Analysis  of  Single  Pings 

Through  the  comparisons  of  single  ping  RL  data  to  RL  averaged  over  the  entire 
event  area,  various  features  can  be  seen  that  are  consistent  in  their  geographic  location. 
The  resulting  backscatter  features  from  Ping  107,  shown  in  Figure  19,  exhibit  RL 
characteristics  that  are  common  in  all  pings  in  the  data  set  and  are  discussed  in  detail 
below.  Figure  19  is  a  geographic  contour  plot  displaying  deviation  of  beam  RL  from  the 
event-wide  mean  RL  for  each  of  the  28  azimuthal  beams  of  ping  107.  It  is  anticipated 
that  regions  where  the  deviations  (from  the  average  RL  curve)  are  positive  correspond  to 
seafloor  areas  exhibiting  more  reflective  characteristics  than  the  area-wide  mean,  while 
regions  with  a  negative  deviation  are  less  reflective.  For  example,  Figure  20  (top)  shows 
the  RL  and  smoothed  RL  curve  of  ping  107,  beam  21 .  At  approximately  0.7  s  a  region  of 
increased  (>5  dB)  RL  is  observed  in  both  the  raw  and  smoothed  curves.  Figure  21 
depicts  this  region  of  ping  107,  beam  21,  enlarged  to  show  the  area  of  heightened  RL 
(2550  m  east,  800  m  north).  As  a  result,  the  deviation  of  beam  RL  from  the  event- wide 
mean  RL  is  referred  to  as  ‘relative-reflectivity’.  The  same  geographic  region  of  increased 
relative-reflectivity  can  be  seen  on  other  contour  plots  using  different  pings.  For 
example,  Ping  97  (Figure  22)  transmitted  about  20  seconds  (~  35  m)  earlier,  and  ping 
117,  (Figure  23)  transmitted  about  20  seconds  (~  35  m)  later,  both  demonstrate  that  the 
region  of  increased  reflectivity  (increased  RL)  is  consistent  in  its  geographic  location 
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(2550  m  east,  800  m  north).  This  indicates  that  the  high  backscatter  or  reflectivity  is  due 
to  a  geographic  feature,  and  is  not  a  random  occurrence. 

Alternatively,  this  feature  can  be  identified  by  a  comparison  between  different 
pings,  by  overlaying  contour  plots  from  a  set  of  pings  from  the  same  event.  Figure  24 
shows  the  sector  areas  covered  by  such  a  sequence  (pings  59,86,113,140,167,  and 
grouped  as  Set  1),  while  their  contour  plots  (Figure  25)  illustrate  how  the  same  high 
backscatter  or  reflectivity  (discussed  above)  appears  consistently  in  the  same  geographic 
location  in  a  number  of  pings,  despite  the  varying  transmission  range  as  the  ship  moved 
progressively  closer  towards  the  features. 

2.  Analysis  of  Events 

A  thorough  analysis  of  all  the  RL  curve  plots  from  each  event  was  first  conducted 
to  create  a  geographic  representation  of  various  sediment  types  based  on  consistent 
spatial  variations  in  the  RL  curves  within  that  event. 

a)  Overlays  of  Ping  Sequences 

Within  an  event  a  continuous  sequence  of  pings,  equally  spaced  at 
approximately  100  m  (4-9  pings  depending  on  ship’s  speed  during  an  event),  was  used  to 
produce  a  display  of  overlaying  contour  plots  that  represents  the  geographic  area  covered 
by  that  particular  event,  as  previously  illustrated  in  Figures  24  and  25.  For  each  event 
this  procedure  was  repeated  four  times,  each  time  with  a  different  set  of  pings,  to  produce 
five  completely  independent  geographic  representations  of  the  area  covered  by  each 
event.  Figures  26-27  show  two  such  sequences  (pings  60,87,1 14,141,169;  Set  2,  and 
pings  61,88,115,142,170;  Set  3,  respectively)  both  from  Event  1  which  demonstrate  their 
consistency  in  portraying  features  evident  in  all  five  representations.  A  comparison  of 
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the  five  displays  from  each  event  reveals  that  the  same  spatial  RL  anomalies  are 
displayed  consistently  in  the  same  geographic  location  regardless  of  the  sequence  of 
pings  used.  A  similar  degree  of  consistency  was  demonstrated  for  other  spatial  RL 
anomalies  in  the  study  area  using  such  sequences  of  RL  ping  data  from  the  other  events 
(Events  2,  3  and  4).  For  example,  the  large  area  of  heightened  backscatter  spanning 
north-south  across  the  sonar  swath  at  1800-2000  m  east  in  Event  2  (Figure  28),  is  not 
only  consistently  displayed  across  all  the  pings  of  that  event,  but  can  also  be  seen  at  the 
end  of  Event  1  (Figure  27,  left)  and  just  at  the  beginning  of  Event  3  (Figure  29,  right).  In 
Event  4  (Figure  30)  several  spatial  RL  anomalies  consistently  occur  near  the  beginning  of 
the  track  (1300-1500  m  east),  but  RL  variation  from  the  mean  appears  consistently  very 
low  in  the  western  end  of  the  area  (0-1200  m  east).  Figure  30  (left)  also  demonstrates 
spatial  consistency  of  RL  features  between  pings  which  progressively  vary  in  centre 
beam  heading  as  the  sonar  was  trained  to  the  right  near  the  end  of  the  event. 

b)  Composite  Plots 

Figure  31  is  similar  to  those  of  Figures  25-27  but  is  a  composite  plot 
derived  from  combining  the  data  of  all  5  pings  from  that  particular  sequence.  For  each 
defined  geographic  location  covered  by  the  sequence,  the  average  deviation  of  RL  from 
the  event-wide  mean  RL  (relative-reflectivity)  was  calculated  and  then  displayed. 

Because  this  is  an  average  of  relative-reflectivity  for  a  number  of  pings, 
this  plot  will  show  RL  features  that  have  geographic  consistency,  while  eliminating  those 
features  that  are  not  geographically  consistent  but  only  appear  in  single  pings.  As  a 
result,  the  composite  plot  appears  to  lose  some  of  the  detail  present  in  the  single  and 
overlay  plots  (Figures' 2 1-27),  but  the  significant  features  are  still  exposed. 
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In  addition,  because  each  ping  from  the  sequence  was  transmitted  at  a 
progressively  shorter  range  from  each  defined  geographic  location,  the  resulting  relative- 
reflectivity  has  been  calculated  from  a  well  distributed  range  of  grazing  angles;  from 
about  1  to  15°  (for  direct  path  transmissions)  as  discussed  in  Chapter  2  (A.l).  This 
effectively  averages  out  any  effect  of  varying  grazing  angle,  and  the  result,  for  each 
geographic  location  of  the  composite  plot,  is  a  generic  value  of  relative-reflectivity, 
independent  of  grazing  angle  (over  this  range). 

Composite  plots  of  the  other  three  events  (Events  2,  3  and  4),  displayed  in 
Figures  32-34,  also  manifest  the  more  relevant  spatial  RL  features  that  were  previously 
displayed  in  their  appropriate  sequence  overlay  plots  (Figures  28-30).  It  should  be  noted 
that  most  of  the  RL  features  observed  in  these  composite  plots  display  nearly  the  same 
level  of  relative-reflectivity  as  they  did  in  the  single  ping  plots  indicating  that  grazing 
angle  plays  only  a  minor  role  in  the  variation  of  RL  (over  the  range  of  shallow  grazing 

angles  (1  to  15  )  analysed  in  this  study).  This  is  in  agreement  with  established  theory  as 
seen  in  Figure  2. 

3.  Combining  Events 

In  order  to  analyse  the  entire  area  of  study  as  a  whole,  the  RL  statistical  data  from 
the  different  events  were  combined  to  produce  an  ‘area-wide’  mean  RL  curve.  However, 
RL  statistical  data  from  one  event  to  the  next  can  only  be  combined  once  the  difference 
between  the  modes  of  operation  (transmission  power  and  pulse  length)  and  their 
implications  on  the  RL  statistics  is  completely  understood. 

When  combining  events  that  employ  the  same  operating  modes,  such  as  Events  1 

and  3,  all  the  RL  statistical  data  from  both  events  can  be  used  to  simply  calculate  an 
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overall  mean  RL  curve  that  represents  the  average  for  both  areas  covered.  However, 
events  which  employ  different  operating  modes  must  be  treated  differently.  For  example, 
a  decrease  in  transmission  power,  or  source  level  (SL),  would  result  in  a  corresponding 
decrease  in  RL  (Urick,  1983).  This  means  that  a  difference  in  SL,  such  as  that  between 
Event  1  and  2  (5  dB),  has  to  be  taken  into  account  when  calculating  the  combined  mean 
RL  for  all  events.  This  involved  applying  an  adjustment  (adding  5  dB  logarithmically)  to 
the  mean  RL  data  of  Event  2,  which  normalised  the  data  (to  medium  power)  so  that  it  was 
equivalent  to  that  of  events  1  and  3,  before  calculating  an  overall  mean. 

The  non-normalised  event-wide  mean  RL  curves  for  Events  1,  2  and  3  are  shown 
in  Figure  35  together  with  a  fourth  which  is  the  mean  of  all  three  mean  RL  curves 
calculated,  as  discussed,  from  all  the  RL  data  of  events  1,  2  and  3.  Notice  that  they  all 
have  similar  shape  and  rate  of  decay,  an  observation  reflected  by  the  correlation 
coefficients  between  one  another  that  range  from  0.99813  (Events  1  &  2)  to  0.9998 
(Events  2  &  3).  This  similar  shape  and  rate  of  decay  is  a  mean  representation  of  the 
environment  and  is  applicable  independent  of  operating  mode  (transmission  power).  The 
displacement  or  offset  between  the  mean  RL  curve  of  Event  2  and  those  of  the  other 
events  is  primarily  due  to  a  difference  in  transmitting  power.  However,  any  RL  level 
difference  between  events  which  had  the  same  operating  mode  (transmitting  power),  such 
as  Events  1  and  3,  can  only  be  attributed  to  a  difference  in  the  environment  (sea  floor 
characteristics). 

The  RL  statistical  data  of  Event  4  was  not  combined  with  that  of  the  other  events 
as  their  relationship  was  more  complicated.  Event  4  not  only  employed  a  different  SL 
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(220  dB),  but  also  employed  a  different  pulse  length,  2  ms  as  opposed  to  1 1  ms.  The 
influence  of  pulse  length  on  RL  can  be  seen  from  the  following  equation  (Urick,  1983): 

RL  =  SL  —  2TL  +  Ss>v  +  logA,  V 

where  TL  =  transmission  loss,  Ss,v  =  scattering  strength  (surface  and  volume),  and  A,V  = 
the  reverberating  area  or  volume,  respectively.  A  is  specifically  given  by: 

A  =  cx/2  x  Or 

where  c  =  sound  speed,  x  =  pulse  length,  O  =  beam  width,  and  r  =  range.  Thus  a  decrease 
in  pulse  length  (x)  such  as  that  between  Events  3  and  4  could  result  in  a  corresponding 
decrease  in  RL,  despite  the  increase  in  SL.  In  addition,  there  was  a  delay  in  time  between 
these  events  (an  hour  and  a  half),  in  which  the  wind  experienced  a  small  but  significant 
decrease  of  2-3  ms  1  (Table  1).  The  resulting  decrease  in  surface  roughness  almost 
certainly  decreased  the  surface  component  of  reverberation,  thus  compromising  one  of 
the  fundamental  assumptions  (stationary  surface  reverberation  over  the  duration  of  the 

experiment),  discussed  in  Chapter  1,  on  which  the  IT  analysis  of  these  combined  events 
depends. 

This  complication  and  resulting  outcome  were  confirmed  when  a  number  of  pings 
from  Events  3  and  4,  transmitted  over  the  same  geographical  area,  were  compared.  The 
88  pings  of  Event  3,  transmitted  at  medium  power  (215  dB),  were  compared  to  the  60 
pings  of  Event  4  (referred  to  as  Event  4a),  transmitted  at  high  power  (220  dB),  which  all 
covered  the  same  geographical  area.  The  event-wide  mean  RL  curves  for  these  two 
events  are  displayed  together  in  Figure  36  where  their  shape  and  rate  of  decay  appear 
identical,  with  a  correlation  coefficient  of  0.99994.  However  the  mean  RL  curve  for 
Event  4a  is  displaced  approximately  2  dB  below  that  of  Event  3,  contrary  to  what  might 
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be  expected  based  on  SL  considerations  alone!  Although  the  mean  RL  from  Event  4a, 
transmitted  at  a  higher  SL,  would  be  expected  to  be  above  that  of  Event  3,  the  effects  of 
both  reduced  pulse  length  and  reduced  surface  reverberation  clearly  had  an  opposite  and 
overriding  effect  on  the  result. 

Analysis  of  these  event-wide  statistics  confirms  that  RL  statistical  data  for  pings 
from  different  events  can  be  combined  to  produce  results  applicable  for  the  entire  area. 
However,  if  the  resulting  measure  of  deviation  (of  single  ping  beam  RL  from  mean  RL) 
is  to  have  any  validity,  each  ping  analysed  must  use  the  appropriate  event-wide  mean  RL 
curve  for  the  event  from  which  the  ping  in  taken. 

4.  Overlay  and  Composite  Plots  of  Entire  Study  Area 

Figure  37  is  an  overlay  of  the  contour  plots  of  a  sequence  of  pings  (Set  1  from 
Events  1,  2  and  3)  covering  most  of  the  study  area.  The  values  of  relative-reflectivity 
have  been  calculated  from  a  single  normalised  area-wide  mean  RL  curve  calculated,  as 
discussed  above,  from  the  RL  statistics  of  Events  1,  2  and  3.  Theses  contour  plots  appear 
to  combine  well,  in  that  overlapping  plots  from  adjacent  events  appear  to  produce  similar 
levels  of  relative-reflectivity  over  the  same  geographic  areas.  Figure  38  is  a  composite 
plot  of  the  same  sequence,  which  is  produced  as  previously  discussed,  but  again  employs 
the  single  area-wide  mean  RL  curve  (Events  1,  2  and  3)  in  calculation.  This  plot  was 
primarily  used  for  the  following  comparison  to  known  sediment  ground  truth. 

5.  Comparison  to  Existing  Information  on  Sediments 

The  previous  sections  have  demonstrated  that  RL  features  consistently  appear  on 
all  of  the  single  ping,  sequence  overlay,  and  composite  contour  plots  produced.  Two 
composite  contour  plots,  one  from  combining  Events  1,  2  and  3  (Figure  38)  and  the  other 
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from  Event  4  (Figure  34),  cover  the  entire  area  and  are  used  in  the  analysis  below  for 
comparison  to  the  (little)  known  sediment  distribution  of  the  area.  However,  it  must  be 
emphasised  that  the  values  of  relative  reflectivity  from  Figures  34  and  38  cannot  be 
compared  directly  to  each  other,  as  the  values  are  relative  to  separate  area/event  mean  RL 
curves,  as  previously  explained.  This  is  confirmed  by  observing  the  difference  in 
appearance,  and  different  values  of  relative  reflectivity,  of  an  area  observed  in  both 
figures  (34  and  38)  at  1300  m  east  and  800  m  north. 

Despite  the  lack  of  spatial  detail  in  all  of  the  existing  sediment  distribution 
sources,  analysis  of  the  contour  plots  show  clear  similarities  to  the  general  sedimentary 
characterisation  of  the  area.  There  are  three  distinct  geographic  areas  or  domains 
observed  in  Figures  34  and  38,  similar  to  that  described  in  Chapter  2. 

The  first  of  these  distinct  domains  is  the  most  significant  feature  in  the  area 
covered  by  Figures  34  and  38,  namely  the  somewhat  featureless  area  west  of  1200  m  east 
and  south  of  1200  m  north.  This  area  displayed  the  least  variation  in  relative-reflectivity, 
as  illustrated  by  a  dynamic  range  of  less  than  1  to  2  dB.  This  implies  a  relatively 
homogeneous  surface/subsurface  sediment  spatial  distribution,  and  is  most  likely  the 
result  of  a  fairly  uniform  surface  sediment  cover  of  mainly  sand,  just  as  the  information 
from  Figure  9  suggests.  The  relatively  smooth  bottom  that  the  bathymetry  implies  in 
Figure  1 1  also  supports  this  conclusion. 

The  second  area  is  the  region  that  extends  from  1200  to  3200  m  east,  although  is 
interrupted  by  the  third  area,  a  band  between  1800  m  to  2100  m  east.  This  second  area 
has  a  high  variability  in  relative-reflectivity  (dynamic  range  as  much  as  10  dB),  and  is 
characterised  by  areas  of  fairly  low  relative-reflectivity  (0.5  to  2  dB  below  the  mean) 
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interspersed  by  irregular  patches  of  quite  high  relative-reflectivity  (2  to  5  dB  above  the 
mean).  This  region  appears  anything  but  homogeneous  and  implies  a  bottom  of  different 
sediment  types.  The  patches  of  high  relative  reflectivity,  such  as  that  around  2400-2600 
m  east  and  600-900  m  north,  and  the  group  around  2700-3100  m  east  and  400-700  m 
north,  all  appear  to  be  quite  localised  sharp  features  of  increased  reflectivity,  some  of 
which  appear  inter-connected.  It  is  suggested  that  these  highly  reflective  zones  could  be 
uncovered  (or  thinly  covered)  rock,  rocky  patches,  or  gravel/pebble  fields,  whereas  the 
less  reflective  areas  may  be  due  to  depressions  or  basins  filled  with  finer  grained 
sediments.  The  few  small  areas  displaying  minimum  relative  reflectivity  (-3  to  -4  dB), 
for  example  1600  m  east  and  900  m  north,  are  quite  likely  to  be  flat  areas  of  fine 
sediment  or  mud.  This  all  closely  follows  the  bottom  type  descriptions  of  this  area, 
especially  those  from  the  MIW  pilot  (Chapter  2),  and  is  indicative  of  the  glacial  origin 
( mixed  moraine)  of  sediments  in  this  area. 

The  third  area  is  by  far  the  most  prominent  feature  in  terms  of  relative  reflectivity. 
The  200-300  m  wide  band  previously  referred  to,  transects  the  entire  swath  of  coverage 
north-south  between  1800  m  to  2100  m  east,  and  stands  out  as  an  extremely  reflective 
feature.  The  reflectivity  and  size  of  such  a  large  feature  suggest  that  it  is  either  an 
exposed  rock  outcrop  or  (ancient)  reef,  or  more  likely  a  particularly  large  linearly 
oriented  field  of  mixed  rock  and  gravel  (moraine).  Although  no  such  feature  is 
specifically  described  in  the  sources  of  sediment  distribution,  the  general  descriptions  of 
mixed  sediment  with  glacial  origin  certainly  make  such  a  feature  a  credible  possibility. 
Interestingly  enough,  this  is  somewhat  supported  by  the  fact  that  the  location  of  this 
feature  is  approximately  coincidental  with  the  area  of  greatest  variation  in  bathymetry 
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(Figure  11),  thus  confirming  that  this  band  has  a  particularly  erratic  or  rough  bottom. 
The  AN/SQQ-32  penetrates  sediment  up  to  approximately  1.5  acoustic  wavelengths 

(Null,  1996),  or  up  to  6  cm.  Thus  partially  buried  rocks  or  outcrops  could  cause  some  of 
the  high  reflectivity. 

In  summary,  the  resulting  geographic  contour  plots  of  relative  reflectivity  not  only 
agreed  with  the  broad  picture  painted  by  the  various  existing  sources  of  information  on 
sediment  distribution,  but  also  provide  far  more  detailed  spatial  representation  of 
sediment  distribution.  This  served  to  highlight  the  gross  inadequacies,  particularly  in 

spatial  resolution,  of  the  existing  information  available  on  bottom  sediment  distribution 
and  geoacoustics. 

B.  DISCUSSION 

1.  Choice  of  Sample  Interval 

All  of  the  analyses  described  thus  far  were  carried  out  with  RL  data  smoothed 
with  a  sampling  interval  of  1000  data  points.  Although  this  interval  is  equivalent  to  a 
range  spacing  of  40  m,  this  does  not  correspond  to  the  actual  spatial  resolution,  as  the 
latter  is  more  complicated  parameter  which  is  also  dependent  on  the  overlap  (75%)!  The 
interval  of  1 000  data  points  was  by  no  means  meant  to  reflect  the  optimum  choice,  but,  as 
mentioned  in  Chapter  2  (B),  was  a  compromise  between  maintaining  sufficient  detail 
(spatial  resolution)  in  the  time  series  data  and  keeping  the  data  sizes  manageable  within 
the  capabilities  of  the  PC  equipment  available  to  this  study.  The  optimum  choice  of 
interval  would  actually  require  further  investigation,  and  also  depend  upon  the  required 
application  and  the  data  management  available.  However,  in  partial  response  to  the 
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question  of  optimum  interval  size  one  single  ping  data  set  was  analysed  using  a  range  of 
different  intervals. 

As  the  sampling  interval  is  directly  related  to  the  spatial  resolution  (and  overlap), 
the  particular  ping  used  for  this  brief  investigation  was  selected  because  it  included  a 
target  of  opportunity  of  known  size.  Ping  21 1  of  Event  2  was  selected  because  a  Mark  36 
mine  (a  cylinder  2.2  m  long  and  0.4  m  in  diameter),  indicated  as  mine  #13  in  Figures  6 
and  7,  had  been  detected  in  this  particular  ping  during  previous  analysis  (Wilson  et  al., 
1996).  It  can  be  seen  in  Figure  39  at  a  range  of  300  m,  predominantly  manifested  in 
beam  1. 

A  series  of  contour  plots  of  ping  21 1  were  produced  and  then  examined  using  a 
progressively  smaller  interval  each  time.  The  intervals  employed  included  2000,  1800, 
1600,  1400,  1200  ,1000,  900,  800,  700,  600,  500,  400,  300,  200,  100,  50,  and  26,  which 
corresponded  to  spatial  intervals  of  80  m  to  1  m.  Examination  of  the  resulting  contour 
plots  revealed  that  even  at  an  interval  as  large  as  2000  (corresponding  to  80  m  range 
interval),  the  broad  picture  of  relative  reflectivity  was  still  very  much  evident.  As  the 
interval  was  decreased,  the  resulting  contour  plots  progressively  revealed  smaller  features 
and  more  detail,  much  as  expected.  For  example,  in  Figure  40  the  contour  plot  of  ping 
21 1  produced  at  an  interval  of  1200  (-50  m)  was  the  first  plot  of  the  series  to  reveal  an 
area  of  high  relative  reflectivity  (-2  dB)  at  2100  east,  500  north.  This  particular  feature 
of  heightened  reflectivity  is  almost  certainly  associated  with  the  quite  significant 
bathymetric  feature,  possibly  a  pinnacle  or  rock  outcrop  of  some  kind,  located  in 
approximately  the  same  position  observed  in  Figure  11.  Mine  13  is  first  observed,  at 


33 


2400  m  east,  470  m  north,  in  the  contour  plot  of  ping  21 1  produced  with  an  interval  of 
200  (8m)  (Figure  41). 

Figure  42  shows  the  contour  plot  of  ping  21 1  produced  using  an  interval  of  26  (1 
m),  the  smallest  interval  used.  This  figure  demonstrates  the  detail  and  high  level  of 
spatial  variation  that  can  be  achieved  using  a  fine  scale  interval,  a  capability  that  would 
prove  most  useful  in  archiving  clutter  densities  and  NOMBO  fixes. 

2.  Analysis  of  RL  Curve  Data  versus  Analysis  of  RL  Slope  Data 

The  same  method  of  analysis  used  on  the  RL  curve  data  was  also  carried  out  on 
the  RL  slope  data,  following  the  procedure  outlined  by  Schalm  (1999)  for  use  on  SQS- 
53C  data.  However  complications  in  removal  of  the  TVG  created  an  (artificial)  artefact 
within  the  RL  slope  data.  This  artefact  corrupted  the  results,  as  illustrated  in  the  contour 
plot  of  Figure  43,  (although  the  most  significant  geographically  consistent  features  were 
still  revealed).  This  corruption  only  occurred  in  the  analysis  of  RL  slope  data,  and  was 
due  to  the  way  in  which  the  TVG  specifically  affected  the  slope  calculations.  The  RL 
curve  data  analysis  was  completely  free  of  such  corruption  as  its  method  involves 
calculations  that  actually  cancel  out  the  artefact.  As  a  result  the  analysis  of  RL  slope  data 
was  put  aside  in  favour  of  the  more  successful  analysis  of  RL  curve  data.  Although  the 
TVG  artefact  was  the  prime  motivation  in  focusing  effort  on  the  analysis  of  RL  curve,  it 
was  firmly  believed  that  this  method  actually  provided  superior  results  anyway,  even  if 
the  TVG  problem  had  not  been  present.  This  however  did  prevent  any  direct  comparison 
of  the  two  methods  (RL  curve  and  RL  slope). 

Comparisons  between  the  results  of  Schalm’s  AN/SQS-53C  work  (1999)  and  this 
study  are  complicated  due  to  the  difference  in  sensors,  applications  and  the  environments 
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involved,  and  thus  one  can  only  make  qualitative  comment.  Notwithstanding  significant 
differences  in  parameters  including  frequency,  pulselength  and  sector  coverage,  one  of 
the  fundamental  differences  between  the  two  applications  is  that  of  spatial  scale;  the 
SQQ-32  extends  a  70°  sector  of  coverage  out  to  a  range  of  approximately  1  km,  while 
SQS-53C  extends  a  120°or  360°  sector  of  coverage  out  to  a  range  of  approximately  19 
km.  Clearly  the  scale  of  area  covered  is  vastly  different.  Further  investigation  involving 
direct  comparison  of  the  two  methods  (RL  curve  and  RL  slope)  applied  to  one  or  both 
applications  (MCM  and  ASW)  is  needed,  but  it  is  expected  that  the  RL  curve  analysis 
will  be  superior  to  the  RL  slope  analyses  for  both  MCM  and  ASW  applications. 

3.  Implications  of  Horizontal  Orientation 

As  the  data  used  was  opportunistic  and  not  gathered  as  a  dedication  effort  for  this 
stud)',  there  were  some  limits  to  the  way  in  which  it  could  have  been  exploited.  One  very 
desirable  aspect  of  this  study  was  not  covered  because  of  this  aspect.  Virtually  the  entire 
area  studied  was  acoustically  surveyed  from  the  same  direction  (transversing  east  to 
west ».  This  prevented  comparison  of  results  using  this  IT  method  over  the  same  piece  of 
bottom  but  from  different  transmission  directions.  It  might  be  expected  that  certain 
types  of  bottom,  particularly  those  with  highly  oriented  features,  such  as  sand  ripples  or 
waves,  would  provide  different  results  depending  on  the  azimuthal  orientation  of  the 
sensor.  To  investigate  this  aspect  would  require  with  tracks  arranged  in  such  a  way  as  to 
obtain  transmissions  over  the  same  geographic  location  from  different  directions,  before 
any  conclusions  could  be  made. 
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4.  Tactical  Implications  of  the  Plots  of  Bottom  Relative  Reflectivity 

The  high  variability  of  bottom  sediment  type  and  clutter  density  in  shallow  water 
presents  a  strong  tactical  significance  because  it  can  unexpectedly,  and  often  times 
adversely,  affect  the  results  of  an  acoustic  search,  especially  when  executed  under  the 
assumption  that  the  bottom  sediment  is  homogenous.  Figure  44(A)  is  indicative  of  a 
typical  MCM  search  plan  that  was  executed  assuming  that  detection  ranges  were  constant 
(regardless  of  sediment  type).  This  "mowing  the  lawn"  technique  can  result  in  lower  or 
higher  probability  of  detections  (PDs)  than  anticipated.  This  is  due  to  regions  where 
detection  ranges  are  degraded  or  improved  because  more  or  less  reflective  areas  are 
encountered  (along  the  search  pattern),  respectively.  Implications  are  two  fold:  If  the  PD 
along  the  search  track  is  lower  than  predicted  the  resulting  search  swath  is  narrower  than 
the  line/track  spacing  planned  for.  This  potentially  results  in  gaps,  known  as  “holidays”, 
between  the  search  lanes,  giving  over-optimistic  confidence  levels  and  inaccurately  low 
risk  estimates.  Conversely,  if  the  PD  is  higher  than  predicted,  the  resulting  search  swath 
is  wider  than  the  line/track  spacing  planned  for,  giving  a  level  of  redundancy,  and 
meaning  that  the  planned  search  takes  longer  than  necessary.  The  overall  implications  to 
a  larger  naval  mission,  highly  dependent  on  the  timely  execution  of  such  a  (MCM) 
search,  are  that  it  proceeds  at  an  unnecessarily  slow  rate,  so  much  so  that  it  may  be 
delayed  or  cancelled,  or  worse,  that  it  proceeds  at  under-estimated  levels  of  risk,  and  thus 
(unknown)  inadequate  levels  of  safety! 

Prior  to  an  acoustic  search  a  rapid  environmental  update  could  be  conducted 
which  includes  an  IT  survey  of  relative  reflectivity,  from  which  a  sediment  map,  similar 
to  the  one  in  Figure  38,  could  be  used  to  modify  or  dictate  the  search  plan.  Such  a 
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survey,  run  at  the  appropriate  spatial  resolution  (sample  interval),  as  discussed  in  Chapter 
2.B.1,  could  also  include  assessments  of  clutter  density  and  even  NOMBO  fixes.  The 
search  plan  could  then  divide  the  area  into  sub-divisions  differentiated  by  PD,  with  the 
line/track  spacing  and  orientation  within  each  sub-division  dictated  by  the  PD  of  that 
particular  sub-division  (Figure  44B).  This  selective  approach,  "searching  by  sediments”, 
would  optimise  the  efficiency  of  such  a  search,  and  areas  of  particularly  high  reflectivity 
(very  low  PD)  or  high  clutter  density,  which  might  be  considered  too  time  intensive  to 
search,  could  be  simply  cordoned  off  and  discarded  in  the  overall  search. 

In  addition,  or  alternatively,  an  IT  analysis  ran  during  the  actual  search  would 
provide  a  very  useful  means  of  real  time  performance  assessment  and  monitoring. 

This  technology  could  be  further  exploited  by  populating  databases,  such  as  the 
Route  Survey  Data  Base  (RSDB),  with  the  actual  processed  IT  data,  either  as  composite 
plots  of  relative  reflectivity  or  in  actual  data  form  such  that  their  plots  can  be  regenerated. 
In  either  case  this  information  would  then  be  available  to  aid  planning  as  discussed 
above,  but  ahead  of  the  arrival  of  fleet  assets.  Such  databases  would  also  facilitate  the 
means  of  comparing  how  bottom  reflectivity  and  clutter  densities  may  have  changed 
between  such  surveys,  a  particularly  useful  attribute  in  respect  to  clutter  densities  and 
NOMBO  fixes.  This  information,  available  to  populate  such  databases,  is  gathered,  and 
currently  wasted,  whenever  a  platform  sensor  is  operating.  This  IT  would  ensure  that  all 
such  data  and  information  is  used  to  its  potential  whenever  possible.  Subsequent 
processing  and  database  population  could  be  made  an  automatic  process. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


USS  Avenger  (MCM-1)  conducted  an  experimental  sonar  trial  in  Rhode  Island 
Sound  on  19  February  1993,  which  included  collection  of  reverberation  level  data  from 
the  ship’s  AN/SQQ-32  sonar.  The  reverberation  level  data  was  used  to  investigate  the 
possibility  of  employing  inversion  techniques  to  infer  bottom  geoacoustic  properties. 

The  hypothesis  presented  in  this  thesis  was  based  on  the  assumption  that  the 
magnitude  of  backscattered  energy  (and  its  rate  of  spatial  decay)  was  directly  related  to 
the  acoustic  reflectivity  of  the  seabed.  By  examining  the  character  of  RL  curves  from 
various  sedimentary  provinces,  one  should  be  able  to  identify  and  associate  various 
reverberant  areas  with  unique  sediment  or  geoacoustic  types.  This  study  examined  the 
statistics  of  both  the  RL  curve  and  its  slope,  to  test  this  hypothesis. 

A.  CONCLUSIONS 

Of  the  two  statistical  parameters  the  RL  curve  was  found  to  be  most  clearly 
aligned  with  the  degree  of  bottom  reflectivity.  Over  less  reflective  (mud/silt)  bottoms  a 
decrease  in  RL  deviation  from  the  area-wide  mean  RL  curve  (relative  reflectivity)  was 
noted  while  acoustic  interactions  with  a  more  reflective  bottom  (gravel,  pebbles  and 
rocks)  cause  an  increase  in  relative  reflectivity.  Based  on  this  finding,  a  technique  was 
dc\ eloped  using  the  deviation  of  the  RL  for  an  individual  ping/beam  from  an  area-wide 
average  RL  to  generate  geographic  maps  illustrating  bottom  acoustic  characteristics. 

Through  the  comparisons  of  different  plots  displaying  the  same  geographically 
located  features  of  reflectivity,  it  was  noted  that  grazing  angle  plays  only  a  minor  role  in 
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the  variation  of  RL  (over  the  range  of  shallow  grazing  angles  (1°  to  15°)  analysed  in  this 
study),  as  in  agreement  with  established  theory. 

Resulting  geographic  maps  of  relative  reflectivity  not  only  agreed  with  the 
existing  descriptions  of  sediment  distribution,  but  also  provided  far  more  detail  in  spatial 
representation  of  the  bottom  geoacoustic  distribution.  This  served  to  highlight  the  gross 
inadequacies,  particularly  in  spatial  resolution,  of  the  existing  information  available  on 
bottom  sediment  distribution  and  geoacoustics. 

In  addition,  geographic  maps  of  relative  reflectivity  can  be  produced,  when  using 
an  appropriately  small  sample  interval,  at  sufficient  spatial  resolution  to  expose  clutter 
density  information  (MILECs,  MELCOs,  NOMBOs  etc.). 

The  resulting  geographic  maps  of  relative  reflectivity  can  provide  the  following: 

1  A  means  of  real  time  performance  assessment  and  monitoring  that  enable 
more  accurate  PDs  and  line/track  spacing,  and  the  ability  to  revise  and 
modify  search  strategy. 

2  An  invaluable  aid  to  planning  a  search  strategy,  optimising  resources 
based  on  “searching  by  sediments”,  either  from  in-situ  IT  surveying  or 
from  archived  IT  derived  data. 

3  A  means  of  comparing  clutter  densities  (MILECs,  MILCOs,  NOMBOs 
etc.)  from  different  surveys. 

4  Appropriate  geoacoustic  parameters  for  propagation  models  (including 
full  wave)  to  provide  more  accurate  predictions. 
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B.  RECOMMENDATIONS 


It  is  recommended  that  the  RL  method  of  this  inverse  technique  analysis  for 
generating  sediment  maps  be  used  by  fleet  assets  as  part  of  efforts  to  develop  both  in-situ 
and  archived  databases  of  bottom  properties,  and  that  these  results  be  incorporated  into 
shallow  water  MCM  (and  ASW)  search  planning. 

In  addition,  further  research  into  the  nature  of  RL  IT  statistics  and  its  impact  on 
sonar  system  performance  should  be  undertaken  and  include  the  following: 

1  Investigation  of  different  sample  intervals  in  respect  to  optimum  spatial 
resolutions  for  both  broad  geoacoustic  distribution  and  clutter  densities. 

2  Investigation  of  IT  results  performed  on  the  same  geographical  area  but 
from  different  azmuthal  directions  to  determine  any  directional 
dependence. 

3  Comparison  of  this  IT  using  both  RL  analysis  method  (Howell,  2000)  and 
slope  of  RL  analysis  method  (Schlam,  1999)  in  respect  to  optimising  both 
ASW  and/or  MCM  applications. 

4  Multi-sensor  data  IT  fusion  algorithms  from  various  sonars  operating  over 
the  same  bottom  area,  but  at  vastly  different  frequencies,  should  be 
developed.  For  US  sonars  the  AN/SQS-53C,  AN/SQQ-32,  and  a  high 
frequency  side  scan  sonar  (SSS)  all  contain  different  information  about  the 
bottom  sediments,  and  their  IT  data  must  be  fused  to  optimise  the  3D 
description  of  the  bottom.  This  equally  applies  to  UK  sonars  including 
Sonar  Types  2050  and  2093. 


41 


It  is  emphasised  that  high  quality  and  accurate  ground  truth  information  is  essential  in 
validating  future  IT  analysis  results. 
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APPENDIX  A.  FIGURES 
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Figure  2.  Bottom  Sb  as  a  function  of  grazing  angle  illustrating  (from  Urick,  1983). 
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Figure  3.  Bottom  backscattering  strength  as  a  function  of  normalised 
grainsize  at  10°  (from  Chotiros,  1992) 
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Figure  4. 


Map  of  Rhode  Island  Sound  displaying  area  of  study. 
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Chart  of  Rhode  Island  Sound  displaying  the  MIW  exercise  area 
(after  NOAA,  1998). 
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Figure  6.  Ship’s  tracks  in  Rhode  Island  Sound  MIW  exercise  area. 
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UTM  (meters  east  of  71 .23  deg  West) 


Figure  8.  The  sediment  distribution  off  the  New  England  continental  shelf  is 


comprised  mostly  of  sand  (after  Milliman  et  al.,  1972). 
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Figure  9. 


Sediment  distribution  in  Rhode  Island  Sound 
(after  NAVOCEANO  MWP,  2000). 
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Figure  10.  Close-up  of  navigation  chart  displaying  ‘Nature  of  Bottom’  descriptors 
showing  sand  throughout  (after  NOAA,  1998). 
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High  Resolution  Bathymetry  of  Area  of  Study  (bottom  contours  in  meters) 
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SQQ32  ARL  Pings  from  Tracks  1  &  2  in  Rhode  Island  Sound  MIW  Exercise  Area 


Figure  12.  Ship  s  ping  distribution  detailing  different  operating  modes 


UTM  (meters  east  of  71.23  deg  West) 


Figure  13.  AN/SQQ32  Time  Varying  Gain 
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SSQ32  ARL  Track  3  Ping  40  at  12:46:15  Beam  14 


SSQ32  ARL  Event  1  Ping  107  at  1 1 :03:44  RL  at  Interval^ OOOpts  &  Overiap=75% 
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Contour  plot  of  Ping  107’s  RL  deviation  from  event- wide  (Event  1)  mean 

RL. 
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SSQ32  ARL  Event  1  Ping  107  at  11:03:44  RL  at  lnterval=1000pts  &.  Overlap=75% 


(gp)  Id  ue3W  UJOJJ  uojiejAaQ 


(IIPOM  Bap  /£'  if  jo  upou  sja;aiu)  ^lfl 


Figure  21.  Region  of  pingl07  beam  21  (enlarged)  to  show  the  RL  curve  details. 
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Figure  22.  Contour  plot  of  Ping  97’s  RL  deviation  from  event- wide  mean  RL. 
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Figure  23.  Contour  plot  of  Pingl  17’ s  RL  deviation  from  event- wide  mean  RL. 
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SSQ32  ARL  Event  1  Set  2  Deviation  of  RL  from  Event  Mean  RL  in  each  Beam  at  I nt e rva 1= 1 000 pt s  &  Overlap=75% 


Figure  26.  Overlaying  contour  plots  from  sequence  of  pings  (Set  2)  from  Event  1. 
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Figure  27.  Overlaying  contour  plots  from  sequence  of  pings  (Set  3)  from  Event  1. 
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SSQ32  ARL  Event  2  Set  1  Deviation  of  RL  from  Event  Mean  RL  in  each  Beam  at  Intervals  OOOpts  &.  Overlap=75% 


(»P/(0P)P)  3(loIS  UB8W  LU0Ji  uo!ie|A8a 


(UPON  Bap  i£  It  P  Wou  sjajaui)  iaiid 


Figure  28.  Overlaying  contour  plots  from  sequence  of  pings  (Set  1)  from  Event  2. 
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Figure  29.  Overlaying  contour  plots  from  sequence  of  pings  (Set  1)  from  Event  3. 
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Figure  30.  Overlaying  contour  plots  from  sequence  of  pings  (Set  6)  from  Event  4. 
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Figure  31.  Composite  plot  of  sequence  of  pings  (Set  1)  from  Event  1. 
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Figure  32.  Composite  plot  of  sequence  of  pings  (Set  1)  from  Event  2. 
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Figure  33.  Composite  plot  of  sequence  of  pings  (Set  1)  from  Event  3. 
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Figure  34.  Composite  plot  of  sequence  of  pings  (Set  6)  from  Event  4. 
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Correlation  between  Events  1  &2  =  0.99813 


Time  (sec) 


SSQ32  ARL  Event  Mean  RLs  at  Interva  1=1000  pts  &.  Overlap=75°/c 


Composite  Relative  Reflectivity  of  SSQ32  Set  51  Deviation  of  RL  from  Event  Mean  RL  at  Intervals  OOOpts  &  Overiap=75% 
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Figure  38.  Composite  plot  of  sequence  of  pings  from  Events  1, 2  &  3. 
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Figure  39  Mine  13  displayed  in  existing  sonar  display 

(from  Wilson  et  al,  1996). 
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Figure  40  Contour  plot  of  Ping  21  l’s  relative  reflectivity  at  interval  of  1200. 
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Figure  4 1  Contour  plot  of  Ping  2 1 1’  s  relative  reflectivity  at  interval  of  200. 

84 


1400  1600  1800  2000  2200  2400  2600 

UTM  (meters  east  of  71. 23  deg  West) 


1200 


(HUOM  Bap  2E  I.fr  P  Mliou  sjaiaui)  iA|in 


Figure  42  Contour  plot  of  Ping  2 1  l’s  relative  reflectivity  at  interval  of  26. 
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Figure  43  Contour  plot  of  Ping  21  l’s  deviation  of  slope  of  RL  from  event¬ 
wide  (Event  1)  mean  slope  of  RL. 
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